Dynamical scaling in Ising and vector spin glasses 
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We have studied numerically the dynamics of spin glasses with Ising and XY symmetry (gauge 
glass) in space dimensions 2, 3, and 4. The nonequilibrium spin-glass susceptibility Xne(£w,T) and 
the nonequilibrium energy per spin e nc (t v , T) of samples of large size Lb are measured as a function 
of anneal time t w after a quench to temperatures T. The two observables are compared to the 
equilibrium spin-glass susceptibility Xcq(L,T) and the equilibrium energy e eq (£, T), respectively, 
measured as functions of temperature T and system size L for a range of system sizes. For any time 
and temperature a nonequilibrium time-dependent length scale L*(t v ,T) can be defined by writing 
Xne(tw, T) = Xcq(L* , T) (or the equivalent expression for the energy). Our analysis shows that for all 
systems studied, an "effective dynamical critical exponent" parametrization L* (t w , T) = J 4(T)i 1//z( - T ' 
fits the data well at each temperature within the whole temperature range studied, which extends 
from well above the critical temperature T c to near T = for dimension 2 or to well below T c for 
the other space dimensions studied. In addition, the data suggest that the dynamical exponent z 
varies smoothly when crossing the transition temperature. 

PACS numbers: 75.50.Lk, 75.40.Mg, 05.50.+q 
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I. INTRODUCTION 

The dynamics of laboratory spin glasses manifests a 
number of fascinating phenomena, linked to the fact that 
below the glass temperature T c the systems never achieve 
true thermodynamic equilibrium. 1,2 It has gradually be- 
come clear that a slow increase of the coherence length 
(also known as "dynamic correlation length" ) with time 
plays an important role in the memory and rejuvenation 
effects seen experimentally in spin glasses under various 
cooling and heating protocolsi&^iSiLS 

Numerical work can also bring light to bear on the 
question. Finite-size scaling theory states that the de- 
pendence of the equilibrium spin-glass susceptibility Xcq 
on the system size L at the critical temperature T c is 
given by 

L 2- n 



(T = T c ), 



(1) 



where r/ is the "anomalous dimension" static scaling ex- 
ponent for the correlation function of the system. In ad- 
dition, at T c , the equilibrium autocorrelation relaxation 
time increases with sample size L as 



-(L) ~ U- 



(2) 



where z c is the dynamical critical exponent, as conven- 
tionally defined for the standard single-spin Glauber up- 
date dynamics. HuseSi remarked that the critical anneal- 
time dependence of the nonequilibrium spin-glass suscep- 
tibility for large samples after a quench to T c is 

Xne(tw) ~ t%-«V Z < , (3) 

where £ w is the time after quench and z c is again the 
equilibrium dynamical critical exponent. Equation J3J is 
strictly equivalent to the definition of an effective time- 
dependent length scale via 

L*{t w ) = A c t l J z ° , (4) 



with A c an appropriate prefactor. The fact that the 
nonequilibrium scaling should depend only on the equi- 
librium dynamical critical exponent z c is important and 
nontrivial. In the first case the system is quenched from 
infinite temperature; thus its effective temperature is 
changing with time all through the annealing process. In 
the second case z c represents the dynamical scaling for 
thermal fluctuations within the set of configurations at 
thermodynamic equilibrium. A rigorous theoretical justi- 
fication for the nonequilibrium approach exactly at criti- 
cality has been given by Jannsen et alM> Nonequilibrium 
dynamics has been studied numerically in considerable 
detail for many regular magnetic systemsfii principally 
because such data can give accurate information on the 
critical behavior (see, for instance, Ref . IT^) . 

Spin-glass critical nonequilibrium behavior has already 
been studied numerically in a number of Ising spin glasses 
(ISG's) and the gauge glass (GG)i 9 i 13 i 14 i 15 i 16 Relaxation 
becomes very slow below T c in glassy systems; the phrase 
"time is length" has been coined for the link between co- 
herence length £(i w ) and anneal time i w r^ and the fact 
has been underlined that the physically relevant length 
scales involved in the dynamics of spin glasses below T c 
are short, even for experimental times which are always 
very long compared to microscopic time scales. This is 
because the values of z c in spin glasses are intrinsically 
high, as the mean field z c at the ISG upper critical di- 
mension d = 6 is already equal to 4, and z c (d) increases 
to yet higher values at lower dimensions. 

In numerical simulations of ISG's below T c the time 
dependence of the dynamical correlation length £(t w ) 
has been estimatedi£*i2*i2i2£ by measuring the time- 
dependent correlation function explicitly and parametriz- 
ing using an appropriate assumption for the form of the 
function. Early correlation length data were analyzed 
using the phenomenological assumption that a dynam- 
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ical scaling relationship of the critical functional form, 
Eq. Q, continues to hold for £(t w ) at temperatures 
below T c , with a temperature-dependent effective dy- 
namical exponent z(T)m&il£ It has been suggested that 

z(T) cx 1/roa 

Alternatively a "dynamic droplet scaling" has been 
proposed, where below T c an excitation barrier increases 
algebraically with correlation lengthiSii 2 ^ For ISG's in 
three and four dimensions the consequences of an analy- 
sis based on this approach have been discussed in detail 
in Refs. IT?! TTflL and Isjj. This form of parametrization 
has been widely employed in analyses of experimental 
data 3 i 6 i 23 i 24 i 25 although it should be noted that length 
scales are never measured directly in experiments. 

In the paramagnetic state (T > T c ), relaxation is 
fast compared with the time scales of most experiments; 
therefore measurements of relaxation rates are very dif- 
ficult. No numerical studies of dynamical length scales 
seem to have been undertaken either in the regime above 
the spin-glass freezing temperature in systems having a 
finite T c , except for the two-dimensional (2D) Edwards- 
Anderson (EA) ISG for which T c = 0, where careful stud- 
ies have been undertakenii2iS& 

In this work we present results from simulations of 
the EA ISG with Gaussian-distributed interactions, as 
well as the GG in dimensions 2, 3, and 4. We choose 
these models as they are paradigmatic representatives of 
spin glasses with Ising as well as with vector spin sym- 
metry. We define a temperature-dependent length scale 
L*(t w ,T) and find that the equation having the same 
functional form as Eq. J3J valid at criticality, 

BL*(t w ,T)=A(T)tl/^ , (5) 

with a temperature-dependent effective exponent z(T) 
and a weakly temperature-dependent prefactor A(T), 
gives an excellent parametrization of the data in each 
system, not only at temperatures below T c (confirming 
the conclusions of Refs. Iisl andll^) . but also above T c . In 
addition, we find a disagreement with "droplet dynamic 
scaling"i2i2£ for the Ising spin glasses, which becomes ap- 
parent for temperatures below ~ T c /2. For temperatures 
above T c neither the droplet approach nor the replica 
symmetry breaking (RSB) approach appears to give any 
predictions as to a dynamic scaling. For both pictures, 
"barriers" in the energy landscape disappear above T c , 
therefore relaxation becoming trivial. 

In Sec. [n] we define and discuss a time-dependent 
length scale which will be needed to rescale the data. In 
Sees. IIIII and II VI we present data for the two-, three-, and 
four-dimensional gauge glass and Ising spin glass, respec- 
tively. After a brief summary (Sec. 01 , droplet dynamic 
scaling is discussed in Sec. IVII followed by concluding 
remarks in Sec. IVIII 



II. SPATIOTEMPORAL SCALING 

The relationship between time and length is defined 
and discussed below using the internal energy and the 
spin-glass susceptibility. In general, the internal energy 
for a spin Hamiltonian is given by 

e=^[(W)]av (6) 

Here N = L d represents the number of spins described 
by a Hamiltonian TL on a hypercubic lattice of linear size 
L and (• ■ •) represents a thermal average, whereas [■ ■ -] a v 
corresponds to a disorder average. 

The spin-glass susceptibility (related to the nonlinear 
susceptibility measured experimentally) can be expressed 
as 

X = N[{\q\ 2 )U . (7) 

Here q represents the Edwards- Anderson order parameter 
which in the case of the gauge glass is given by 

1 N 

? =-5;eipW-^)], (8) 
J=l 

with <fij representing the phases of the XY spins, and a 
and (3 are two copies of the system with the same disor- 
der. For the Ising spin glass q is given by 

1 N 

where Sj = ±1 represent Ising spins. 

A. Definition of a length scale 

For any fixed T, as the system size L is increased or 
the anneal time i w lengthed, the SG equilibrium and 
nonequilibrium susceptibilities Xeq(L,T) and Xnc(tw,T) 
grow while the internal equilibrium and nonequilib- 
rium energies per spin, e cq (L,T) and e ne (t w ,T), re- 
spectively, drop towards the infinite-size equilibrium 
value e eq (oo,T). If T is above the ordering tem- 
perature, Xeq(L,T) and Xne(t w ,T) will both saturate 
at a temperature-dependent limiting value x cq (oo,T); 
e eq (L,T) and e ne (t w ,T) will always saturate at eoo(T) 
for all T. 

Quite generally, if at any temperature T the equilib- 
rium SG susceptibility as a function of sample size L is 
Xeq(L,T) and the nonequilibrium SG susceptibility after 
an anneal time t w following a quench to temperature T is 
Xne(tw, T) for a large sample of size Lb, then an infinite- 
sample-size time-dependent length scale L*(t w ) can be 
rigorously defined by writing 

Xeq [L*(t w ),T]= X ™(t w ,T) , (10) 
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as long as L*(t w ,T) <C Lb- This length scale defini- 
tion applies to any temperature and any anneal time, 
and the length scale L* (t w , T) can be estimated numer- 
ically to high precision if data of sufficient statistical 
accuracy are available. A practical limitation to preci- 
sion under some conditions is the need to intrapolate 
between the Xeq(L, T) sequence of values for integer L 
to provide a continuous function with which to compare 
the time-dependent data. In the paramagnetic regime, 
as Xcq(L,T) and \ne(tw,T) both saturate at the equi- 
librium infinite-size Xeq(oo,T) value, L*(t w ,T) can only 
be determined up to some temperature-dependent finite 
time. 

An analogous equation can be derived for the energy 
per spin, 

e eq [L*(t w ),T]=e ne (t w ,T), (11) 

providing an independent estimate for L*(t w ,T). Later, 
we show that similar estimates of L*(t w , T) are obtained 
from a completely independent analysis of the data for 
the susceptibility x an d energy e. 

Note that L*(t w ,T) is not the time-dependent coher- 
ence length £(t w , T)>il but it is closely related to it. Dur- 
ing the anneal each individual spin is surrounded by a 
growing cohort of spins in equilibrium correlation with 
it at the temperature T, up to a time-dependent cutoff 
length. This correlated volume can be considered equiv- 
alent to a box of linear size L*(t w ,T). Because £(t w ,T) 
is defined as a correlation length and L*(t w ,T) from the 
box size, we should expect L*(t w ,T) ~ 2£(i w ,T). Here, 
£(t w ,T) has been estimated from directly measured cor- 
relation functions, although a fully rigorous definition is 
not easy to give because of nontrivial prefactorsiiiiS*i2iSli 
For the 3D EA ISG where the comparison can be made 
directly, the equivalence between L*(t w ,T) and 2£(t w , T) 
is indeed found to be correct. It should be borne in mind 
that because of the disorder, the length values represent 
averages over different samples and the effective length 
scale L* (i w , T) may be slightly different depending over 
which observable the mean is being taken. 

B. Comparison between effective length scales of 
different observables 

There is a straightforward way to test the assump- 
tion that the observables Xnc(tw,T) and e nc (t w ,T) are 
controlled by a single length scale. At each T, the 
equilibrium energy e cq (L,T) can be written as a func- 
tion of the equilibrium susceptibility Xeq(L,T) with L as 
an implicit parameter: e eq (L,T) = f e [x oq (L, T)}. Sup- 
pose that Eqs. 1)10(1 and 1(11(1 hold, with one and the 
same L*(t w ,T) for both observables. Then the nonequi- 
librium energies e nG (i w ,T) are given by the identical 
function f e [Xne(t w , T)} of the nonequilibrium suscepti- 
bility Xne(tw,T), with t w an implicit parameter (see 
Fig.nj. A plot of nonequilibrium [e nc (f w , T) — e cq (oo, T)] 
against Xnc{tw,T) data is superimposed on a plot of 



[e eq (L,T) - e oq (oo,T)] against Xcq{L,T) data [by defi- 
nition the equilibrium and nonequilibrium energy mea- 
surements tend to the identical e eq (oo,T) which we esti- 
mate by extrapolation] . These figures are simple displays 
of raw data, and no fitting procedure whatsoever is in- 
volved; at this stage no assumption is made as to the 
functional form of L*(i w , T). 

In what follows we present data for the gauge glass, as 
well as the Ising spin glass. 

III. GAUGE GLASS 

The gauge glass is a canonical vector spin glass (see, 
for instance, Refs. where XY spins on a (hy- 

per)cubic lattice in d dimensions of size L interact 
through the Hamiltonian 

H = -jJ2 cos (<Pi-4 , j-A ij ), (12) 

(i,3) 

the sum ranging over nearest neighbors. The angles 
<pi represent the orientations of the XY spins, and the 
Aij are quenched random variables uniformly distributed 
between [0, 27r] with the constraint that Aij = —Aji 
(here J = 1). Periodic boundary conditions are applied. 
The GG ordering temperatures have been shown to be 
T c = 0, 0.46(1), and 0.89(1) in dimensions 2, 3, and 4, 
respectivelyi 16 i 28 i 30 i 31 i 32 

For XY spin systems there is a choice to be made 
in the allowed single-spin acceptance angle for individ- 
ual updating steps. To optimize the updating procedure 
at low temperatures, the limiting angle is often chosen 
to be less than 2ir for an XY spinS and linearly de- 
pendent on temperature. The numerical prefactor for 
the temperature-dependent window is chosen so that the 
acceptance ratios for the local Monte Carlo updates is 
~ 0.4. As far as the final equilibrium parameters are 
concerned, this choice plays no role. However, for the 
nonequilibrium simulations it is essential to use the full 
2ir acceptance angle window to obtain physically signifi- 
cant results as otherwise the limited angle introduces an 
artificial temperature variation in the relaxation. 



A. Two dimensions 

The GG in space dimension 2 has a zero-temperature 
ordering transitioni 29 i 30 i 31 i 32 i 34 i 35 i 36 i 37 Therefore all the 
measurements discussed in this section necessarily refer 
to the paramagnetic state. Dimension d = 2 presents the 
advantage that systems can be equilibrated up to large 
L, so in at least part of the temperature range compar- 
isons can be made between Xne(^w, T) and Xcq(L, T) over 
a wide range of system sizesi^ Details of the simulations 
are presented in the Appendix, Table ITTI for equilibrium 
and Table lyTTTl for nonequilibrium measurements, respec- 
tively. 
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As the systems are always paramagnetic at finite T, for 
the whole temperature range Xeq(L) must finally saturate 
at the thermodynamic infinite-size limit and a plot of 
\og[xcq(L,T)] against log(L) is always curved (although 
the curvature is weak at small sizes and low T). At 
each temperature a comparison between the nonequilib- 
rium data for the susceptibility Xnc(t w ,T) and energy 
Zne(tw,T) and the equilibrium susceptibility Xeq(L,T) 
and energy e cq (L,T) can be made. Figures show 
examples at two different temperatures T = 0.173 and 
T = 0.513. In each case a plot is first made of 

e p (i w ) = [ene(tw, T) - e oq (oo, T)) (13) 
against Xnc(t w ,T) and of 



e p (L, T) = [e cq (L, T) - e oq (oo, T)} 



(14) 



against Xcq{L,T). The excellent superpositions show 
that to good accuracy the effective length scale L*(t w ,T) 
for the two observables is the same throughout the an- 
neal at each temperature. In both equations e eq (oo, T) is 
estimated by extrapolation. The evaluation of L*(t w ,T) 
from the energy data is not sensitive to e eq (oo, T) as long 
as exactly the same value is used for the equilibrium and 
nonequilibrium data. 

To obtain L*(t w ,T) explicitly, the Xae(t w ,T) data are 
scaled onto the equilibrium Xeq{L, T) data, translating to 
x(L*,T) by assuming the space-time relation in Eq. (J5J 
and adjusting A(T) and z(T) to obtain optimum scaling. 
Scaling plots are shown for the same two temperatures. 
The condition L*(t w ,T) <C L\> is well satisfied for all 
the data with the present ranges of sizes and maximum 
anneal time. At the higher temperatures studied (as in 
Fig.Hfor T = 0.513), Xcq{L,T) and Xnc{t w ,T) saturate 
and the range of points from which the fit is usable is 
restricted to lengths and times before the onset of satu- 
ration. For T = 0.513 accurate scaling can still be ob- 
tained, but this condition leads to a practical upper limit 
on the temperatures over which A(T) and z(T) can be 
estimated. Note that these two temperatures are chosen 
as typical examples of the behavior in the two regimes 
in which the data sets do not or do arrive at satura- 
tion, respectively, within the available ranges of L and 
t w . In other space dimensions and for the Ising systems 
the equivalent plots always take on one or other of the 
behaviors depending on the temperature. 

The size-dependent equilibrium energy data 
e p (L,T) — [e e q(L,T) — e(oo,T)] and the time- 
dependent nonequilibrium energy data e p (t,T) = 
[e ne (t w ,T) — e(oo,T)] translated to e p (L*,T) can also 
be plotted together assuming exactly the same scaling 
as for the susceptibility data as discussed above. To the 
present accuracy, for these temperatures the z(T) values 
obtained from the energy data are indistinguishable 
from the values estimated from the susceptibility scaling 
[although the effective prefactors A(T) are slightly 
different]. The errors in each z(T) data point are 
subjective estimates obtained by varying the parameters 
around the optimal values. 
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FIG. 1: (Color online) Energy difference between the en- 
ergy per spin for a large sample (L = 128) at time f w af- 
ter quench and the infinite time energy, Eq. (11311 . plotted 
against the large-sample spin-glass susceptibility after quench 
Xne(tw) (triangles) for the two-dimensional GG at T = 0.173. 
On the same plot the equilibrium energy difference for size 
L [Eq. 11411 ] is plotted against the equilibrium susceptibility 
Xeq(L) for the same size (squares). These are raw data, and 
no fitting is involved in these plots. 
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FIG. 2: (Color online) Finite-size equilibrium susceptibilities 
Xcq(i) vs L (squares) for the two-dimensional GG at T = 
0.173. In addition, the susceptibilities Xne(tw) at time i w after 
a quench are plotted against L*{t„) — A(t„)tl/ Z ' T ' with the 
parameters A and z, both temperature dependent, adjusted 
for optimal overall overlap with the equilibrium data points 
(triangles). In this figure L*(t„) = 1.66tw . 
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FIG. 3: (Color online) Finite-size equilibrium energy differ- 
ences e p (L) [Eq. 1141 j vs L (squares) for the two-dimensional 
GG at T = 0.173. On the same plot the energy differences 
e P (t w ) at time i w [Eq. I13H 1 after a quench are plotted against 
L*(t w ) = A(t v ,)tlf z< - J " > with identical parameters A(T) and 
z(T) as used in Fig. [5] (triangles). In this figure e p = e+ 1.502. 
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FIG. 5: (Color online) Two-dimensional GG equilibrium 
(squares) and nonequilibrium (triangles) susceptibility data 
plotted as in Fig. |2] for T = 0.513 with appropriate adjust- 
ment for optimum A(T) and z(T) for this temperature. In 
this figure L*(t v ) = 1.346£w 08 . 
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FIG. 4: (Color online) Two-dimensional GG equilibrium 
(squares) and nonequilibrium (triangles) susceptibility and 
energy data for T = 0.513, as in Fig. Q Note the nontriv- 
ial functional form. 




FIG. 6: (Color online) Two-dimensional GG equilibrium 
(squares) and nonequilibrium (triangles) energy data are plot- 
ted as in Fig. for T = 0.513 with the same parameters A(T) 
and z(T) used in Fig. In this figure e p — e + 1.263. 



The estimates for z(T) against 1/T are shown in Fig- d 
The data can be parametrized using 

z(T) = 2.18(9) + 0.41(2)/T, (15) 

with x 2 — 0.60 for the fit. This implies a diverging dy- 
namical exponent z(T) as T approaches zero and z(T) 
tending to near 2 at infinite temperature. Note that z = 2 
for a random walk. 

We conclude from this section that for the GG in di- 



mension 2 which orders at zero temperature, a time- 
dependent length scale L*(t w ,T) can be measured over 
a wide range of temperatures in the paramagnetic state. 
This length scale obeys the effective dynamical critical 
exponent scaling relationship, Eq. (JSJ. The sets of values 
of z(T) determined from the nonequilibrium spin-glass 
susceptibility Xnc{twiT) and the nonequilibrium energy 
per spin e ne (t w ,T) are the same to within the present 
precision. 
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FIG. 7: (Color online) Effective dynamical exponent z(T) as a 
function of 1/T for the two-dimensional GG. The dashed line 
corresponds to a fit of the form z(T) = 2.18(9) + 0.41(2) /T 
and is meant as a guide to the eye. 



Three dimensions 




FIG. 8: (Color online) Effective dynamical exponent z(T) as a 
function of 1/T for the three-dimensional GG. The dashed line 
corresponds to a fit of the form z(T) = 1.95(8) + 1.17(4) /T. 
The shaded region denotes T c — 0.46(1). In this and the 
following figures, the width of the shaded area represents the 
error bar on the estimate of the critical temperature. 



The GG in dimension 3 has an ordering temperature 
T c = 0.46(l)ii£i2& The analysis protocol used is essen- 
tially the same as in dimension 2 (Sec. lIIl'X|) : parameters 
of the simulation are listed in the Appendix, Table ITTT1 for 
equilibrium and Table IVIIII for nonequilibrium measure- 
ments, respectively. Below T c and in the L and t w ranges 
that we have studied, the equilibrium susceptibility in- 
creases as Xeq(L, T) oc L X ( T ^ and the nonequilibrium sus- 
ceptibility as Xnc{tw,T) oc tw with T-dependent x(T) 
and y(T). This algebraic behavior facilitates the analysis 
because the log-log plots of \ are all straight lines. 

As in the two-dimensional case, at each temperature 
an effective dynamical exponent z(T) and prefactor A(T) 
can be defined from the scaling of the equilibrium and 
nonequilibrium susceptibility data using Eq. JSJ. As in 
two dimensions the energy data can be scaled satisfac- 
torily using the same z(T) obtained at each T from the 
analysis of the susceptibility data. Below T c the prefac- 
tors A(T) are slightly different for both the energy, as 
well as the scaling of the susceptibility (not shown). 

In three dimensions z(T) also increases ~ 1/T and has 
an intercept z(oc) w 2 — i.e., 

z{T) = 1.95(8) + 1.17(4)/T, (16) 

with x 2 — 1-06 for the fit. 

The present estimate for z c at T c = 0.46(1) is z c = 
4.5(1), in agreement with previous estimates^ 

It is important to note that z(T) traverses T c smoothly 
with no apparent anomaly; see Fig. [H] This implies that 
in terms of the evolution of length scales with time, the 
dynamics above and below the ordering temperature fol- 
low the same pattern although in the final equilibrium 



configuration all the spins are ordered below T c and while 
they are only correlated over a finite length scale above 
T c . In this system the critical behavior is not exceptional 
as far as the length scale dynamics is concerned. 

C. Four dimensions 

The GG in four dimensions has an ordering transition 
at T « 0.89(1).— We perform a similar analysis as done 
in the two- and three-dimensional case. Parameters of 
the simulation are listed in the Appendix, Tables llVl and 
I Villi A very similar pattern of behavior is observed as in 
the GG with lower space dimensions, with z(T) increas- 
ing linearly with inverse temperature — i.e., 

z{T) = 1.75(13) +2.4(1)/T. (17) 

Here x 2 — 0.3 for the fit. 

As in dimension 3 there is no sign of any change of 
behavior in the region of T c ; see Fig. El At T c w 0.89 we 
obtain z c = 4.5(1), in agreement with Ref. 1 1 (i 

IV. ISING SPIN GLASS 

The Hamiltonian of the Edwards- Anderson Ising spin 
glass 3 — is given by 

H = — JjjSjSj, (18) 

where the sum is over nearest-neighbor pairs of sites on a 
hypercubic lattice in d dimensions, the Si are Ising spins 
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FIG. 9: (Color online) Effective dynamical exponent z(T) as a 
function of 1/T for the four-dimensional GG. The dashed line 
corresponds to a fit of the form z(T) = 1.75(13) + 2.4(1)/T. 
The shaded region denotes T c = 0.89(1). 




FIG. 10: (Color online) The effective dynamical exponent 
z(T) as a function of 1/T for the two-dimensional ISG. The 
line is a guide to the eye. 



taking values ±1, and the are Gaussian distributed 
with zero mean and standard deviation unity. Simula- 
tions are done using periodic boundary conditions. Pa- 
rameters of the simulation are listed in the Appendix. 



A. Two dimensions 

As in the case of the two-dimensional GG, the EA ISG 
with Gaussian interactions in dimension 2 orders only at 
zero temperature ^SJLMJ^JUk thus, all the data con- 
cern the paramagnetic regime. Details of the simulations 
are summarized in the Appendix, Tables IVl and HXl The 
data show that an analysis according to Eq. (jSJ) provides 
an excellent parametrization of the growth of correla- 
tions. The results are in good agreement with direct 
measurements of the correlation functions 

The energy per spin data, e(i w ,T), can be 
parametrized consistently in terms of the same effective 
set of dynamical exponents z(T) as used for the analysis 
of the x(T) data. 

The temperature dependence of z(T) is much stronger 
than for the GG systems, and z(T) deviates some- 
what from a linear variation with inverse temperature; 
cf. Fig.UHl For 1/T < 5, the z{T) data can be approx- 
imately parametrized by z(T) ~ 3.9/T. (It is possible 
that the curve could bend at low 1/T to a nonzero in- 
tercept) The present analysis is consistent with those of 
Refs. Hi and 111 who obtain z(T) ~ 4.3/T using a differ- 
ent numerical technique. In Fig. [23 the dashed line is a 
quadratic fit in i = 1/T and should serve as a guide to 
the eye. 



B. Three dimensions 

There is general consensus that the freezing temper- 
ature of the three-dimensional E A ISG with Gaussian 
interactions is T c = 0.94(2) (Refs. I46l47l48l49|) and that 
the dynamical critical exponent is z c ~ 6.45(10). The 
equilibrium SG susceptibility Xoq(T,T) and energy per 
spin, e cq (T,T), are measured at temperatures between 
T = 0.1 and T = 3.0 for several intermediate sample 
sizes, see Table IVII Nonequilibrium measurements are 
made for Lb = 24; see Table llXl For temperatures up to 
T c and for the range of sizes used, to a good approxima- 
tion (as in the GG below T c ) 

Xcq (L,T) = C(T)L*^ , (19) 

with C(T) a temperature-dependent prefactor and x(T) 
an effective exponent which becomes equal to the true 
static critical exponent 2 — r\ at T c . x(T) can never exceed 
d; for this system it reaches values very close to 3 as 
T tends to zero. C(T) remains very close to 1 for the 
whole temperature range. We can note that by definition 
Xeq(L = 1) = 1 for all T. In the low-temperature range, 
for the equilibrium susceptibility x oq (T,T) we make the 
assumption that log[x C q(T, T)] extrapolates linearly with 
log(T) to zero at L = 1. (L = 2 measurements have not 
been used, as at this particular size there are intrinsic 
"wrap-around" problems associated with the definition 
of the interactions.) 

The effective dynamical critical exponent scaling 
parametrization analysis using Eq. @ is satisfactory over 
the whole temperature range covered. Below T c , the 
present zjT) values, Fig. ^2 are consistent with those 
of Refs. and \lQ but more accurate. The prefactors 
A(T) vary only slightly with T, and A(T) ~ 2B(T) where 
B(T) are the prefactors for the coherence length £(t w , T) 
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FIG. 11: (Color online) The effective dynamical exponent 
z(T) as a function of l/T for the three-dimensional ISG. The 
line corresponds to a fit of the form z(T) = 6.40(15)/T. The 
squares correspond to data for z(T) as estimated from the en- 
ergy, whereas the triangles represent data for z(T) calculated 
using the spin-glass susceptibility. The inset shows the pref- 
actor A of the scaling relation in Eq. @ as estimated for the 
energy and susceptibility. The data show little temperature 
dependence. 



estimated by Refs. Ha and fl9t This agreement confirms 
the conjecture made above that as a general rule the cor- 
relation length £(t w ) can be taken as equal to L*(t w )/2 
to a good approximation. For temperatures above T c , a 
scaling of the effective exponent form remains very sat- 
isfactory; see Fig. ^2 

As found by Refs.lTslandllli z(T) varies approximately 
linearly with l/T. The present data (which are more ac- 
curate than those of the previous work) , including points 
in the paramagnetic region up to about 2.5T C , are con- 
sistent with z(T) ~ 6.40(15)/T ( x 2 ^ 0.85). The points 
for z(T) vary smoothly and continuously through T c , as 
in the three- and four-dimensional GG. 

Again, one can carry out a scaling plot for the nonequi- 
librium energy e ne (t w ,T) in the same way as for the sus- 
ceptibility. The effective z(T) values estimated from the 
energy scaling are consistent with the values from the 
susceptibility scaling, but the prefactors A(T) become 
slightly different in the lower temperature range. This 
confirms that one single anneal-time-dependent length 
growth law controls both susceptibility and energy dur- 
ing the anneal, which seems a more satisfactory form of 
analysis than, for instance, that given in Sec. VI A of 
Ref.M 

It would be of interest to carry out further measure- 
ments so as to obtain significantly more information at 
low and moderate temperatures, but this would require 
nonequilibrium simulations extending to much longer an- 
nealing times. 




FIG. 12: (Color online) The effective dynamical exponent 
z(T) as a function of l/T for the four-dimensional ISG. The 
line corresponds to a fit of the form z(T) — 9.15(20)/T. 



C. Four dimensions 

For the EA ISG with Gaussian couplings in four di- 
mensions it is known that T c = 1.78(2); 15 i 18 i 47 with a 
dynamical critical exponent z c = 4.9(2)M^ The present 
data are analyzed in just the same way as for the three- 
dimensional ISG. There is excellent agreement between 
estimates for z(T) between susceptibility and energy 
data. 

In four dimensions z(T) varies approximately linearly 
with inverse temperature — i.e., z(t) ~ 9.15(20)/T with 
X 2 — 2.8; see Fig. The estimate for z at T c is z c = 
5.1(2), in agreement with previous workii^ Once more, 
z(T) evolves smoothly as T passes through T c . 



V. SUMMARY 

The observed behavior of the dynamical exponent z(T) 
for the six systems studied is summarized in Table [I] It 
can be seen that for the GG systems the data in each di- 
mension can be parametrized in the form z(T) ~ a + b/T 
with the constants a ~ 2 and b increasing with dimension 
d. For the EA ISG systems, z(T) ~ b/T—i.e., a = 0, 
with b values again increasing strongly with space di- 
mension d; in each dimension the value of b is higher 
than that of the GG by a factor of 5-10. For each of 
the ISG's the data could be compatible with a nonzero 
infinite-temperature intercept if there is some curvature 
in z(l/T) at temperatures higher than those we have 
studied. In two dimensions the data potentially suggest 
such a behavior. 
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TABLE I: Temperature dependence of the dynamical ex- 
ponent z(T) — a + b/T for the different models studied at 
different space dimensions d. GG refers to the gauge glass, 
ISG to the Ising spin glass with Gaussian interactions. For 
the Ising spin glass in two dimensions the quoted value of b 
corresponds to a linear fit which is a poor approximation to 
the behavior seen in this system (marked with an asterisk). 



System da b 

GG 2 2.18(9) 0.41(2) 

GG 3 1.95(8) 1.17(4) 

GG 4 1.75(13) 2.4(1) 

ISG* 2 3.9(1) 

ISG 3 6.40(15) 

ISG 4 9.15(20) 



VI. DROPLET DYNAMICS APPROACH 

We have obtained very satisfactory scaling of the Ising 
spin glass and gauge glass dynamics using the effec- 
tive exponent parametrization, Eq. 0. We now discuss 
an alternative analysis based on the dynamical droplet 
approach, which distinguishes critical behavior near T c 
from activated behavior with a barrier exponent ip fur- 
ther below T c ml Berthier and Bouchaudii^S (BB) carried 
out four-spin correlation function measurements on the 
three-dimensional EA ISG similar to those of Refs.lLSIand 
liH but analyzed their data on the basis of the droplet 
dynamics formula^ 

A similar analysis of Heisenberg spin-glass dynamics 
(see Fig. 11 of Ref. Ifjlj) shows qualitatively different be- 
havior from that of the systems reported here. We have 
no data on Heisenberg systems although the present tech- 
nique could potentially be applied there as well. 

The BB droplet dynamics leads to a growth law with 
time corresponding to a coherence length £(t) 

t w [£(t,T)]^T*[£(t,T)]^exp[T(T)£(t,Tf/T} , (20) 

where z c is the dynamical critical exponent, ip is the bar- 
rier exponent, T(T) = T (l - T/T c )'^, and r* is a pref- 
actor. BB scale the data using z c = 7, v — 1.65, tp = 1, 
and To = 5.5 and assume r* ~ 2 but to find good overall 
fits they needed to choose temperature-dependent values 
for t*. 

Yoshino et alM made a closely related analysis of mea- 
surements on the four-dimensional ± J ISG. They intro- 
duce a "crossover length" 

L {T) = T l ^{l-T/T c )- u (21) 

and a "crossover time" 

r (T) = [Lo(T)]*« (22) 

together with a general scaling law for the coherence 
length, 

£(t w )/L (T) = f e [t w /r Q (T)} , (23) 



which relates the coherence length £(t w ) to the elapsed 
time t w . Here fe represents an unknown scaling function. 
The scaling curve depends on the critical parameters (T c , 
z c , and v) which can be taken as known from measure- 
ments at T c , and on the barrier exponent tp. It can be 
readily shown that the BB equation, Eq. Q2U[I. is an in- 
version of the scaling equation of Yoshino et al., Eq. H23|) . 
with an explicit functional form which can be written as 

t w /r (T) = r*[£(t w )/L a (T)P exp{T [^ w )/L (T)^} . 

(24) 

We now analyze the present data for the EA ISG in 
three dimensions in terms of the BB /Yoshino droplet dy- 
namics approach; 52 In Fig. ^] we plot the present data 
in three dimensions for L*(t w ,T) as a function of i w 
at different T in the scaling form proposed by Yoshino 
et ai.)Sa i.e., we plot L*(t w ,T)/2L (T) against t w /T (T) 
with L (T) = T 1 /^! - T/T c )~" and r = [L (T)] Z < [the 
factor 1/2 corresponds to the translation from L*(t w ,T) 
to £(ty,,T)] (see inset of Fig. Illfl . For this scaling we 
adopt for the three-dimensional EA system the param- 
eters proposed by BB: z c = 7, v = 1.65, and ip = 1. 
According to Yoshino et al., the entire data set scaled 
this way should lie on a single (but unspecified) scaling 
curve. The extra BB scaling parameters define one spe- 
cific scaling curve within the same plot. On the "Yoshino 
plot" we therefore draw the full BB scaling plot, Eq. J23J, 
using the two remaining fit parameters from the BB pa- 
rameter set: To = 5.5 and r* = 2. 

In the temperature region from T c down to about T c /2 
the overall Yoshino scaling and the agreement between 
the present data and the BB scaling curve is acceptable. 
This is the same temperature range as covered by the 
BB simulations. However, in the temperature range be- 
low T c /2 where the activated droplet dynamics should 
be valid because little is affected by critical dynamics, 
the scaling breaks down. The curves for different T are 
not superimposed, and the deviations from the BB curve 
correspond to many orders of magnitude along the time 
axis. 

An analysis of the EA ISG data in four dimensions 
leads to the same conclusion. The measurements of 
Yoshino et al~ were made on the ±J model while 
the present measurements correspond to the model with 
Gaussian disorder, but if we adopt the critical parameters 
T c = 1.8, z c = 4.5, and v = 0.93, together with ip = 4.5 
from Yoshino et ai.^ we obtain an overall scaling plot 
for the present EA ISG results in four dimensions. Once 
again the agreement between the scaling prediction (a 
unique scaling curve) and the data set is acceptable for 
temperatures T c /2 < T < T c , but it breaks down at lower 
temperatures (not shown). 

Thus for these two canonical ISG systems the standard 
droplet dynamics scaling below T c does not give a satis- 
factory global account of the data; by suitably adjusting 
the various fit parameters the droplet scaling relation, 
Eq. H23fl . can be made to conform reasonably well to the 
data over limited ranges of temperature and annealing 
time, but at low temperatures this form of scaling pre- 
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FIG. 13: (Color online) Three-dimensional ISG data plotted 
following to the parametrization of Yoshino et al. (Ref. l2fT) . 
According to this droplet dynamics approach all the data 
presented in this way should lie on a single scaling curve. 
The long-dashed line is the particular scaling curve calcu- 
lated using the specific parameters proposed by Berthier and 
Bouchaud in Ref. ^3 f° r this system. For the definitions of 
the parameters Lo and to see the text, Eqs. 12111 and . 
respectively. 

diets a time scale t/ro which diverges much more rapidly 
with increasing length scale L* /Lo than is observed. 

If we attempt to determine ip for the three-dimensional 
GG from scaling our data using this droplet expression, 
we find only poor scaling whatever tp values are assumed 
for a very wide range of values. 

Kisker et al.^S. as well as fitting very successfully their 
EA ISG data in two and three dimensions to Eq. JSJ, 
also compared their data to a different phenomenological 
droplet-dynamics- inspired formula: 

£(t w )-£ = X(T)[ln(t w )] 1 ^ (25) 

and found that this relation was capable of giving a sat- 
isfactory fit to their low-temperature data. However, as 
they point out in Ref. Il9l this formulation has a num- 
ber of defects; in particular, one would not expect the 
droplet formula to be applicable in the paramagnetic re- 
gion. We can also examine the mathematical basis of 
the equation. Suppose that the true time dependence 
is £(t w ) = Ati /z —i.e., \n[£(t w )} = ln(A) + (l/z)\n(t w )- 
with z large. The parametrization, Eq. 125|) . can be re- 
cast as£(t w ) = £ -(l + x) or ln[l(t w )] = Ln(4) + ln(l+a:), 
with x = [Blntfv,)] 1 ^, where B = {\{T)/£ )^ is small. 
Because ln(l + x) = x — x 2 /2 + we have ln[£(t w )] = 
\n(£ ) + [B\n(t w )} 1 ^-[B\n(t w )} 2 ^/2 + --, If > 1, 
on a log[£(t w )} against log(f w ) plot there will be an up- 
bending from the second term compensated by a grow- 
ing negative third term. Trial and error shows that if 
ip ~ 2/3, this parametrization produces to quite good ac- 
curacy a linear dependence of log[£(t w )} on log(i w ) over 



a wide range of t w , typically from t w = 10 to t w = 10 8 , 
successfully mimicking ln[£(t w )] = ln(A) + (1/z) ln(i w ). 
Thus a fit to the parametrization in Eq. I|25[) can be ob- 
tained with a pseudo "barrier exponent" ip even when the 

1 / z 

true behavior is better described by £{t w ) — AtJ , so this 
fit is a mathematical artifact. The value of ip obtained 
from such a fit can be expected to be ~ 2/3 regardless of 
the system studied, which explains why the apparent ip 
estimated in Ref. [l^ both for the three-dimensional ISG 
well below T c and for the paramagnetic two-dimensional 
ISG at low temperatures are close to this value. 

We conclude from this section that the standard 
droplet scaling approach does not give a satisfactory over- 
all account of the numerical Ising spin-glass dynamical 
data in dimensions 3 and 4. 



VII. CONCLUSION 

We have studied the dynamical behavior of ISG's and 
GG's in dimensions 2, 3, and 4 as functions of anneal 
time t w after a quench to a temperature T. We define a 
time-dependent correlation length scale L*(t w ,T) by re- 
lating the time-dependent nonequilibrium susceptibility 
Xne(t w ,T) to the size-dependent equilibrium susceptibil- 
ity Xcq(L,T) and the time-dependent energy e nc (£ w ,T) 
data to the equilibrium energy e cq (L,T) data. In each 
system and at all temperatures (below and above the 
freezing temperature T c , as well as at T c ) a growth 

law for the length scale L*(t w ,T) = A(T)tl/ z(T} gives 
a good parametrization of the data, with an effective 
temperature-dependent dynamical exponent z(T) and a 
prefactor A(T) which is only weakly temperature depen- 
dent. In each case, independent estimates of z(T) from 
the susceptibility and from the energy measurements are 
the same within the precision of the measurements; the 
prefactors A(T) corresponding to the susceptibility and 
to the energy can be marginally different. z(T) diverges 
approximately as 1/T at low temperatures in all the sys- 
tems, and at high temperatures it appears to tend to 
a limiting value consistent with z(oo) ~ 2 for the GG 
systems and to a value which extrapolation suggests is 
close to zero for the ISG systems. The evolution of z(T) 
as the temperature passes through the freezing temper- 
ature T c is smooth for all four systems having a nonzero 
T c . In the Ising spin glasses the dynamic droplet criti- 
cal scaling approachiiiSi is incompatible with the present 
measurements for the region below T c /2. The droplet 
picture gives no predictions whatsoever concerning time- 
length relationships in the paramagnetic region above 
T c . The data show that the standard droplet dynam- 
ics assumptions^ii22i22iS of an effective barrier height in- 
creasing as a power of £ and disappearing at T c are inap- 
propriate for the systems studied. 

Many recent experimental measurements of dynamics 
in spin glasses^S*2iS have been interpreted using the 
droplet dynamics parametrization. It is nontrivial to 
transpose conclusions obtained from numerical data to 
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the experimental regime; the effective time scales for ex- 
periments are vastly greater than for the numerical re- 
sults, and length scales are slightly larger. It must also 
be kept in mind that there appears to be no way in which 
to measure spin-glass coherence lengths directly from ex- 
periment. As the present results invalidate the stan- 
dard droplet dynamics in Ising spin glasses at least on 
time scales attainable numerically, the question should 
be raised as to what physical significance, if any, can 
be ascribed to parameters such as the barrier exponent 
ip which are obtained from scaling analyses of experi- 
mental data based on standard droplet dynamics. As 
has been pointed out in Ref. 0, for one canonical ISG 
system the range of published estimates of -0 from dif- 
ferent experiments extends from 0.03 to i^Q^^^^ 
which also suggests that the standard droplet dynam- 
ics parametrization is inappropriate. Recently, a rigor- 
ous bound for the barrier exponent inferred from two- 
dimensional calculations^ has been given with 0.25 < 
tp < 0.54. Still, it would seem well worth attempting to 
review the experimental data to see if they can be rein- 
terpreted in terms of an effective exponent length scale 
growth scenario compatible with the present numerical 
results. 

The effective dynamical critical exponent scaling sce- 
nario with z(T) oc 1/T which provides a satisfactory 
parametrization for the ISG systems is very similar to 
that observed by Paul et ulMt-i-L for random-bond Ising 
ferromagnets (RBIM's) and diluted ferromagnets in di- 
mension 2. The physical mechanism for the time de- 
pendence of the correlation length in these ferromagnetic 
models is domain coarsening, and the data are inter- 
preted in terms of effective barriers to domain growth 
which increase logarithmically with sizeS It would seem 
very plausible to ascribe the correlation length growth in 
spin glasses to an analogous mechanism. Equation (9) 
of Paul et aJM^L can be written as z(T) = 2 + c/T, c a 
constant, in our terminology, which is precisely what we 
have observed empirically for the GG systems. However, 
it should be noted that at the ordering temperature T c 
for the RBIM the effective dynamical exponent tends to 
the pure system value z ~ 2i£& In the ISG's and GG's 
the dynamical critical exponent z c is always much higher 
than 2 and z(T) continues to decrease regularly as the 
temperature is raised through the paramagnetic regime. 
This is true both above the ordering temperature when 
T c is nonzero or at all temperatures when T c = which 
is the case in dimension 2. 

The standard approach to the dynamics of the growth 
of the coherence length £(t w ) with anneal time after 
quench t w is to assume that there are three principal 
relaxation regimes: paramagnetic at temperatures well 
above T c , critical in the region around T c , and activated 
at temperatures well below T c , each regime having a 
qualitatively different relaxation behavior. The present 
data show that for Ising spin glasses and for the gauge 
glass, a dynamical relationship of the standard critical 
form, £(t w ) = A(T)t\j ' z<yT ^ with a temperature-dependent 



TABLE II: Parameters of the equilibrium simulations for the 
two-dimensional gauge glass. -/V samp is the number of samples, 
iVsweep is the total number of Monte Carlo sweeps for each of 
the 2Nt copies (two replicas per temperature) for a single 
sample, and Nt is the number of temperatures used in the 
parallel tempering method. The lowest temperature used is 
0.13, the highest 1.058. For L = 24 the lowest temperature 
studied is 0.20. 



Nt 

3 10 000 8.0 x 10 4 30 

4 10 400 8.0 x 10 4 30 
6 10 150 8.0 x 10 4 30 
8 8 495 2.0 x 10 5 30 
12 6 890 8.0 x 10 5 30 
16 2 500 2.0 x 10 6 30 
24 2 166 2.0 x 10 6 24 



dynamical exponent z(T), gives a good account of the 
nonequilibrium dynamics at each temperature and not 
only at T c . The data indicate that z(T) varies smoothly 
as a function of temperature when passing through T c . 
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APPENDIX: NUMERICAL DETAILS 

Equilibrium measurements are carried out with sam- 
ples fully thermalized using the exchange Monte Carlo 
(parallel tempering) techniquei^SiSi To ensure that the 
system is equilibrated, we perform a logarithmic data 
binning of all observables (energy and spin-glass suscepti- 
bility) and require that the last three bins logarithmically 
spaced agree within error bars and be independent of the 
number of Monte Carlo sweeps (MCS) iV SWCO p. In the 
case of the Ising spin glass we use the equilibration test 
for short-range spin glasses first introduced in Ref. Irj2l 
Details about the equilibrium simulations for the gauge 
glass are summarized in Tables ILTl 11111 and II VI for d = 2, 
3, and 4, respectively. Details about the equilibrium sim- 
ulations for the Ising spin glass are presented in Tables 
IVirvTl and lvTTl for d = 2, 3, and 4, respectively. For all 
runs we ensure that the parallel tempering Monte Carlo 
moves have acceptance probabilities of at least 30%. 

Parameters used in the nonequilibrium simulations are 
summarized in Tables IVIIII and llXl for the GG, as well as 
the ISG, respectively. 
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TABLE III: Parameters of the equilibrium simulations for 
the three-dimensional gauge glass. The lowest temperature 
simulated is 0.05, the highest 0.947. The different quantities 
are explained in the caption of Table [H] 



L N Eei mp A swocp Nt 

3 10 000 6.0 x 10 3 53 

4 10 000 2.0 x 10 4 53 

5 10 000 6.0 x 10 4 53 

6 5000 2.0 x 10 5 53 
8 2000 1.2 x 10 6 53 



TABLE IV: Parameters of the equilibrium simulations for 
the gauge glass in four dimensions. The lowest temperature 
used is 0.70, the highest 1.345. The different quantities are 
explained in the caption of Table [H] 



L 


-/Vsamp 


-/Vswccp 


Nt 


3 


5000 


2.0 x 10 4 


17 


4 


5000 


8.0 x 10 4 


17 


5 


5000 


4.0 x 10 5 


17 



By convention, values quoted at time t correspond to 
an average taken between MCS's [t + 1] to [2t] following 



an anneal of t MCS's. 



TABLE V: Parameters of the equilibrium simulations for the 
Ising spin glass in two dimensions. The different quantities 
are explained in the caption of Table [H] For L — 128 the min- 
imum temperature used is 0.20. For all other system sizes the 
minimum temperature is 0.05. The maximum temperature 
used is 1.391. Note that for L < 16 standard parallel tem- 
pering Monte Carlo is used, whereas for L > 32 the cluster 
method by Houdayer (R.ef. l63j) is applied. 



L 


-A^i samp 


-iVj sweep 


N T 


3 


10 000 


2.0 x 10 5 


20 


4 


10 000 


2.0 x 10 5 


20 


6 


10 000 


2.0 x 10 5 


20 


8 


10 000 


2.0 x 10 5 


20 


10 


10 000 


1.0 x 10 6 


20 


12 


10 000 


1.0 x 10 6 


20 


16 


10 000 


1.0 x 10 6 


20 


24 


10 000 


1.0 x 10 5 


20 


32 


10 000 


1.0 x 10 5 


20 


64 


1000 


1.0 x 10 6 


40 


128 


250 


1.0 x 10 6 


63 
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